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Abstract. We have solved for the potential flow down- 
stream of the terminal shock of the solar wind in the 
limit of small departures from a spherical shock due to 
a latitudinal ram pressure variation in the supersonic 
solar wind. The solution connects anisotropic stream- 
lines at the shock to uniform streamlines down the he- 
liotail because we use a non-slip boundary condition on 
the heliopause at large radii. The rotational velocity 
about the heliotail in the near-field solution decays as 
the fourth power of distance from the shock. The polar 
divergence of the streamlines will have consequences for 
the previously discussed magnetic pressure ridge that 
may build-up just inside the heliopause. 

Introduction 

We have previously shown the conditions under which 
potential flow occurs beyond the terminal shock of the 
solar wind. Topologically, potential flow represents well 
many of the global properties which are universal to 
all heliospheric stagnation point flow topologies [Swess, 
1990]. The flow is self-consistently incompressible to 
about 15% and irrotational to 10% [Swess and Nerney, 
1990, 1991 (SN1, SN2)]. The latter is true when the 
non-slip condition is used on the heliopause at large 
radii. Because it is likely that reconnection occurs near 
the stagnation point between the heliosheath and inter- 
stellar magnetic fields, we believe that the heliopause is 
best represented as a non-slip interface [5ue,ss and Ner- 
ney, 1993]. This reduces the solar wind Mach number, 
M, down the heliotail compared with free-slip numerical 
models. Then the effect of compressibility in this region 
is reduced by M 2 . Recently, a numerical simulation for 
an external Mach number of M=0.8 has shown that 
the shock is nearly spherical and Sun-centered even in 
that extreme case [, Steinolfson , 1994]. That calculation 
showed that the nose-to-tail density in the heliosheath 
changes by a factor of two; our use of a non-slip he- 
liopause greatly reduces this effect. 

We do not consider the case of external supersonic 
flow because we believe that the interstellar wind is 
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likely to be sub-Alfvenic and information will be able to 
propagate upstream, removing the necessity for a bow 
shock (also see the discussion in Nerney et al . , 1995). 

The merging of the shocked, subsonic solar wind with 
the 25 km/s interstellar wind was calculated in our po- 
tential models by including the effects of a finite radius 
(R s ) terminal shock. This effect entered the solution 
for the streamlines through a parameter e which is a 
measure of the ratio of the interstellar ram pressure to 
that of the pressure in the shocked solar wind just be- 
yond the terminal shock. Vanishingly small values of 
e reproduce Parker’s [1963] model for a weak interstel- 
lar medium with a large turning radius for the solar 
wind. Because of this, the shock radius is small com- 
pared to the transverse dimension of the heliosphere, 
and the shock appears to have a negligible radius. The 
numerical solutions were plotted for variable 6, and the 
distance from the stagnation point to the terminal shock 
was no larger than R s for probable heliospheric param- 
eters. 

It is our purpose now to show the nature of the change 
in these results when a heliocentric latitudinal ram pres- 
sure variation is imposed on the supersonic wind up- 
stream of the shock. This requires a non-spherical shape 
for the shock due to, for instance, high-speed streams 
emanating from the poles of the sun, and changes the 
character of the solutions near the shock. Observations 
of a latitudinal variation in the interplanetary Lyman 
a distribution [ Kumar and Broadfoot , 1979; Lallement 
et al ., 1985] have been interpreted as indicating a 30- 
50% decrease in the solar wind mass flux from solar 
equator to pole during 1973-1977. This would not lead 
to a large equator-to-pole gas pressure variation on the 
shock for our model. The Bernoulli equation analysis in 
SN1 (p.6409) shows that the pressure just downstream 
of the shock only depends on the pressure at large dis- 
tances down the tail on each streamline to about 6%. 
Because we have no velocity shear across the tail at 
large radii, the pressure is the same on all streamlines 
at large r and equal to that in the interstellar medium, 
P ioo . Because we have linearized the solutions in the de- 
parture from sphericity, we are limited to the low end 
of the equator-to pole mass flux decrease. 

The mapping from heliocentric spherical polar coor- 
dinates to heliocentric heliotail coordinates (cylindrical) 
introduces an azimuthal dependence of the streamlines 
about the heliotail. There is also rotational flow about 
the heliotail which decays as the fourth power of dis- 
tance from the shock. A further asymmetry is intro- 
duced because the footpoints for each streamline begin 
at different distances from the sun. 
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Mathematical Formalism 


We solve the incompressible, irrotational (potential 
flow) equation for fluid flow: 

V 2 $ = 0. (1) 


so that the angular variables occur in both the numera- 
tors and denominators of equation (7), unlike the spher- 
ical shock solutions. We have defined e as in SN2: 


The velocity is derived from -V<$ subject to the bound- 
ary conditions: 

e r * v(9, <f>) \ Rs - v 9 {0,<j>) (2a) 

lim v — Vimiz (26) 
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where R s is the distance to the terminal shock, v s is the 
solar wind flow velocity downstream of the shock, <f> is 
the azimuthal angle about the heliotail, and 0 is the po- 
lar angle with respect to the z-axis, the direction of the 
interstellar wind. We impose the following asymptotic 
flow pattern on the supersonic solar wind: 


Now equating equations (3c) and (7) and linearizing 
the various powers of \/V using binomial expansions, 
allows us to solve for the coefficients of the spherical 
harmonics so that 4> may be written as: 
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where 9s is the solar polar angle and 6 < .2, and v s o is 
the solar wind velocity downstream of the shock along 
the solar rotation axis. Rotating into heliotail, helio- 
centric spherical polar coordinates: 

v s — v s q[1 - tf(cos 2 0 -f sin 2 0cos 2 <£)]. (36) 

The rotation of coordinate systems introduces a cos 2 <j> 
dependence on azimuth about the heliotail. Equiva- 
lently, 

V , = W,o[l - jr(4Po + 2f> 2 + P 2 cos 20)] (3c) 

where Pj and P ( m are the Legendre polynomials and the 
associated Legendre polynomials, respectively. 

Now, <I> may be written as: 
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The velocities are derived from the derivatives of <2>. 
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and R 3 q is the value of R s along the solar rotation axis. 
The Am y n are determined by requiring that the velocity 
of the shocked wind approach that of the interstellar 
medium at large r, eq.(2b). As in SN 1 
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In passing, we note, as expected, that there are no order 
e effects in v <p. 

We will now switch to cylindrical coordinates aligned 
with the heliotail and convert the derivatives in (11) to: 


which is derived by setting the ram pressure of the solar 
wind equal to the net pressure in the local interstellar 
medium (LISM). Eq. (6) serves to define V and gives 
the angular variation for the distance to the termination 
shock. 

Taking the radial derivative of equation (4), including 
only terms present in u 5 , and evaluating at R s , we find: 
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For economy of notation we have written r 2 for zu 2 + z 2 . 
The streamlines are defined by 


dz dzv zud<f> 
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Using the nondimensional variables from SN1 
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(note that e 1 / 3 / R s o remains finite as e — ► 0) the stream- 
line equations can be written as: 
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and Fi, F<i are rewritten as 
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When 7] = 0 and dr]/d( is infinite, v z must be zero. 
Setting the left denominator of (15) to zero then deter- 
mines the non-dimensional distance to the stagnation 
point: 

Cc(l - 6F 2 ) + (C 2 ) 3/2 - e(l - 2.45) = 0 (19) 


Numerical Analysis 

Figure 1 shows the streamlines for a 20% deviation 
from a spherical shock both in the meridional plane 
(solid lines; <f) — 90°, which includes the solar rotation 
axis) and in the equatorial plane (dashed lines; <f) = 0°). 
The equatorial streamlines map back to the inner circle, 
a cross-section of the terminal shock, while the merid- 
ional streamlines map back to a cross-section of a pro- 
late spheroid and show the 20% bulge over the poles. 
The positions of the heliopause for the polar and equa- 
torial streamlines are indicated on the right edge of the 
graph, showing that the bulge in the shock over the 
poles pushes the heliopause out in the meridional plane 
compared to the equatorial plane. The streamlines are 
nearly radial on the outside of the shock for <j> = 90°, 
but there is elevated non-radial flow on the shock for 
<f> ~ 0° and e — .125 (an appropriate value for the he- 
liosphere). We choose to superimpose the graphs of the 
streamlines for = 0",90° because these are the two 
angles about the heliotail for which v $ is zero (see eq. 
(lid)). The following plots indicate that the problem is 
inherently three-dimensional for any other value of <f>. 

Figures 2a, b, and c show three views of streamlines 
that begin on the terminal shock at <f> = 60°, followed by 
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Figure 1 . Figure 1 . A plot of streamlines for the equa- 
torial plane [<f> — 0°; dashed lines and circular shock) 
and the meridional plane (<f> ~ 90°; solid lines and oval 
shock for a 20% elongation over the poles). We use 
€ = .125 - an appropriate choice for the heliosphere. 
Arrows on the right-side of the figure (A for <j> = 0°, B 
for 90°) point to the respective heliopauses. Streamlines 
are fully three-dimensional for other <f>. 


fig 2d which superimposes both 4> ~ 30° (grey curves) 
and (j> = 60° (black curves). The effects of v<f> are seen in 
the near-field solutions for the potential flow but decay 
asymptotically as r ~ 4 (see eq. (lid)). 

The potential flow solutions connect the axisymmet- 
ric boundary condition at infinity (all streamlines blend 
into the assumed uniform interstellar wind) with the 
anisotropic boundary condition on the terminal shock 
(high-speed streams at the poles) so that the stream- 
lines must diverge away from the poles more than at 
the equator. Another way to describe this effect is to 
remember that the solution to Laplace’s equation mini- 
mizes the total kinetic energy per gram integrated over 
any volume bounded by streamlines. This requires a 
divergence of the flow away from the poles where the 
kinetic energy per gram on the shock is high. 

Summary and Conclusions 

We have solved a potential flow model for the flow 
down stream from a non-spherical terminal shock. The 
equations were linearized for small departures from 
sphericity and the streamline equations were numeri- 
cally integrated to show the three-dimensional flow ef- 
fects. The anisotropic streamlines on the terminal shock 
diverge away from the polar regions, while in the dis- 
tant heliotail they smoothly take on the characteristics 
of the assumed uniform streamlines in the interstellar 
medium. This would have consequences for the previ- 
ously reported ridge of magnetic flux that may build- 
up just inside of the heliopause [Nerney, Suess, and 
Schmahl, 1993], The divergence of the flow away from 
the poles will initially open magnetic flux tubes, re- 
ducing the magnetic flux in the ridge. We have also 
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Figure 2. Figure 2. Panels (a), (b), and (c) show different views of the streamlines that begin 
on the terminal shock at <j> = 60°, while (d) shows a superposition of <f> — 30° streamlines (grey 
curves) and 4> — 60° streamlines (black curves). Only panels (b) and (c) include streamlines in 
the interstellar medium. In panel (a) the view is looking directly at the sub-stagnation point on 
the terminal shock, while in (b) and (c) the stagnation point is to the right. The streamlines 
diverge from the poles as the anisotropic boundary condition on the shock (high-speed streams 
at the poles) evolve into the uniform flow of the interstellar wind in the downstream heliotail. 


reported the magnetic field line topology and solar cy- 
cle imprint for a kinematic magnetic field in the helio- 
sphere [ Nerney, Suess, and Schmakl , 1995]. A non- 
spherical shock would change the detailed shapes of 
both the magnetic field lines and the solar cycle im- 
print shown in that calculation. However, the largest 
difference shows up in using a free-slip versus a non-slip 
(our condition) boundary condition on the heliopause. 
The consequence of a free-slip boundary condition is to 
increase the flow speed down the heliotail, with a corre- 
sponding decrease in the cross-sectional area of the tail 
and in the thickness of the helioshcath on the upstream 
side of the heliosphere. 

To the extent that velocity shear exists in the he- 
liosheath inside the boundary layer at the heliopause, 
higher speed flow of the shocked solar wind in the in- 
terior of the heliosheath might lead to the importance 
of compressibility, which we have not included. We es- 
timated the effects of compressibility in our model at 
about 15% (SNl), but we cannot be sure of the global 
effect with a non-slip boundary condition on the he- 
liopause until a non-linear model is generated. 
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